Symmetry of boundary conditions of the Dirac equation for electrons in carbon 

nanotubes 
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We consider the effective mass model of spinless electrons in single wall carbon nanotubes that 
is equivalent to the Dirac equation for massless fermions. Within this framework we derive all 
possible energy independent hard wall boundary conditions that are applicable to metallic tubes. 
The boundary conditions are classified in terms of their symmetry properties and we demonstrate 
that the use of different boundary conditions will result in varying degrees of valley degeneracy 
breaking of the single particle energy spectrum. 



INTRODUCTION 

Carbon nanotubes are the subject of intense research, 
motivated by the desire to use their unique physical and 
electronic properties in the development of nanoscale 
electrical devices [1, 2]. The electronic properties of nan- 
otubes follow from the band structure of graphene - a 
two-dimensional (2D) sheet of graphite - which is a semi- 
metal, having a vanishing energy gap at the six corners, 
K-points, of the hexagonal first Brillouin zone. A single- 
wall nanotube may be thought of as a graphene sheet 
rolled up to form a nanometre-diameter cylinder. Period- 
icity around the circumference results in quantized trans- 
verse wavevectors leading to metallic or semiconducting 
behaviour depending on whether the K-point wavevector 
K is an allowed wavevector. 

While the energy spectrum of an infinitely long tube 
will be continuous, a finite tube should possess discrete 
energy levels corresponding to standing waves typical of 
a confined quantum particle. Evidence of discrete levels 
was seen in transport measurements [3, 4] a few years ago, 
followed by the direct observation of sinusoidal standing 
wave patterns by scanning tunneling microscopy [5, 6]. 
The measured wavelength of the standing waves A ~ 
0.75nm, about three times larger than the lattice con- 
stant a w 0.25nm, corresponded to wavevectors near the 
K-point K. More recently Coulomb blockade measure- 
ments on carbon nanotube quantum dots [7, 8, 9] have 
found evidence for fourfold periodicity of the spectra that 
is in agreement with expectations based on spin and K- 
point degeneracy, although the experiments appeared to 
show varying degrees of degeneracy breaking. A number 
of authors [10, 11, 12, 13, 14] have modelled finite-length 
nanotubes in order to describe the atomic scale varia- 
tion of standing waves patterns and the opening of an 
energy gap displaying an oscillating dependence on the 
tube length. Rather than concentrating on one partic- 
ular model of a boundary, we aim to describe all pos- 
sible energy independent hard wall boundary conditions 
for metallic single wall nanotubes. We will classify the 
boundary conditions in terms of their symmetry proper- 
ties and show how different boundary conditions produce 
varying degrees of K-point degeneracy breaking. 



In the scanning tunneling microscopy measurements of 
Ref. 6 an additional slow spatial modulation of the stand- 
ing waves was observed. It was interpreted as being a 
beating envelope function with wavevector q, |q| <C |K|, 
resulting from the interference of left and right moving 
waves with slightly different total wavevectors K ± q. 
Theoretically, the effective mass model [15, 16, 17, 18] 
provides a reliable analytical description of the electronic 
structure near the K point where the total wavevector is 
k = K+q and the dispersion relation is linear E = sv |q|, 
v is the Fermi velocity and s = ±1 for the conduction and 
valence band, respectively. For spinless electrons, the en- 
velope wavefunction \P (q, r) has four components corre- 
sponding to two inequivalent atomic sites in the hexago- 
nal graphite lattice ( "A" and "B" ) and to two inequiva- 
lent K-points in the hexagonal first Brillouin zone. The 
resulting eigenvalue equation for \t is the massless Dirac 
equation, 



Zen.VV = £X': '» = f ',. " ! : (1) 
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where the role of spin ( "pseudo-spin" ) is assumed by the 
relative amplitudes on the A and B atomic sites: a is a 
vector in the (x, y) plane rotated by the chiral angle r/ 
of the tube. Also, v = (V3/2) 07 is the Fermi velocity, 
a is the lattice constant of graphite and 7 is the nearest 
neighbour transfer integral. 

In this paper we consider the effective boundary con- 
ditions for the envelope function ^ in a finite size carbon 
nanotube. Since the effective mass model for \& corre- 
sponds to the Dirac equation, we begin by deriving all 
possible energy independent hard wall boundary condi- 
tions for the Dirac equation. We write them in terms of 
a small number of arbitrary parameters, mixing angles, 
that describe mixing between boundary conditions with 
different discrete symmetries. Then, in order to illustrate 
the meaning of the general boundary conditions, we eval- 
uate the resulting energy level spectra for non-interacting 
electrons in finite-length metallic nanotubes and the cor- 
responding standing wave envelope functions. To antici- 
pate a little, we find that energy independent hard wall 
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boundary conditions for the Dirac equation may be ex- 
pressed in general terms as 



* = M*; 



M 



1; 



{n B .a,M} = 0, (2) 



where M is an Hermitian, unitary 4x4 matrix M 2 = 1 
with the constraint that it anticommutes with the oper- 
ator riB-a, proportional to the component of the current 
operator normal to the interface, riB is the unit vector 
normal to the interface. As explained in the Appendix, 
we find four possible linear combinations of matrices sat- 
isfying these constraints on M, which, assuming ne is 
a vector confined to the (x, y) plane, may be written in 
terms of a small number of arbitrary parameters: 



Mi 
M> 
M 3 
M 4 



cos A (7 n ® ni. a) + sin A (U z <g> n 2 .cr) , (3) 
cosT(!/i.n®7 <7 )+sinT(f2.n®n B .c7), (4) 
cosfi (^2-n <g> ne-cr) + smSl (In ® ni .a) , (5) 
cos6 (v 1 .Il®I a )+ sin6 (n z ®n 2 .c7), (6) 

where the angles A,T,e and f2 are arbitrary, ni and 
n 2 are three-dimensional space-like vectors satisfying the 
contraints ni.ne = n 2 .riB = ni.n 2 = 0, and v\ and u 2 
are two-dimensional (confined to the (x, y) plane) space- 
like vectors satisfying the constraint u^.u^ = 0. Here we 
have adopted a matrix direct product notation to high- 
light the separate K-point space and AB space structure, 
using the notation {a x , a y , a z , I a } for 2 x 2 Pauli matri- 
ces and the unit matrix that operate within a block ('AB 
space') and {H x , n y , H z , In} for 2 x 2 Pauli matrices and 
the unit matrix that operate in K-point space. For ex- 
ample, the operator a may be written as a direct product 
a = H z ® a. Note that the boundary conditions of the 
Hadron bag model[19], in which elementary particles are 
confined by a scalar (mass) term at the boundary, are 
described by M = v^.Yl ® n^.a. Berry and Mondragon 
[20] considered "neutrino billiards" with a two compo- 
nent Dirac equation confined by a term proportional to 
a z , corresponding to M = In ® ni.cr or M = U z ® n 2 .cr 
with cither ni or n 2 lying in the (x, y) plane. 

There are two non-equivalent K-points that we label 
as K and K. The Dirac equation is diagonal in K-point 
space, so that, in the absence of boundary conditions, 
there are two right moving (^^ and \&^) and two left 

moving (*^ } and *^ } ) plane wave solutions near the 
Fermi surface of a metallic tube. The solutions and 
are eigenvectors of the pseudo-spin component along 
the tube axis S.ne with eigenvalue +s, whereas the so- 
lutions and have eigenvalue — s, where s = ±1 
for the conduction and valence band, respectively. Also, 
the solutions and are eigenvectors of pseudo- 
helicity — iS.V/|q| with eigenvalue +s, whereas the so- 
lutions \J>^ and have eigenvalue — s. There arc dif- 
ferent ways of combining the right and left moving waves 
in order to create standing waves. The first possibility 



is that waves at the same K-point combine, namely 
and form a standing wave with helicity eigenvalue 
+s, and and form a standing wave with helic- 
ity eigenvalue — s. Tms situation is realised by the matrix 
Mi, Eq.(3), because it is diagonal in K-point space. A 
second possibility is that waves from opposite K-points 
combine, namely and form a standing wave 

with spin component eigenvalue +s, and \I>^ and \£'^ 
form a standing wave with spin component eigenvalue 
— s. This situation is realised by the matrix M 2 , Eq.(4), 
because it is off-diagonal in K-point space. A third pos- 
sibility is a combination of the previous two, with waves 
scattered back at the boundary into a mixture of both 
of the K-points. This situation is realised by the matri- 
ces M3, Eq.(5), and M4, Eq.(6), because they have both 
diagonal and off-diagonal in K-point space parts. 



EFFECTIVE MASS MODEL 

In the effective mass model of two-dimensional 
graphite [15], the total wavefunction ^ to t is written as 
a linear combination of four components m = {1, 2, 3,4} 
corresponding to two K-points (i = {1,2} and 7r-type 
atomic orbitals <fj(r — Rj) on two non-equivalent atomic 
sites j = {A, B} in the unit cell, 

* tot (r) = (r) - G ro (r) .V + . . .} i> m (r) , 



m—1 



(7) 



where 



« = -^E^-^r-R,.), ( 8 ) 
1 N 

G ro (r) = -^^e^^tr-RjKr-Rj), (9) 



R, 



arc Bloch type functions constructed from the atomic or- 
bitals, Rj is the position of an atom in real space and the 
summation is over the number of unit cells N ^> 1. The 
functions ip m (r) are components of the envelope function 
^ (q, r) . Substituting this expression for ^ tot into the 
Schrodinger equation and integrating with respect to fast 
degrees of freedom that vary on the scale of the unit cell 
leads to the Dirac equation Eq.(2) for the envelope func- 
tion ^ where the K-points are taken as K = (±47r/3a, 0) 
and the components of ^ are written in the order KA, 
KB, KB, KA. The appearance of the chiral angle of the 
tube 77 in the Dirac equation shows that the axes of the 
(x, y) coordinate system have been rotated to be trans- 
verse and parallel to the tube axis. Applying periodic 
boundary conditions to the wavefunction ^tot, Eq.(7), in 
the direction transverse to the nanotube axis produces a 
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condition for the envelope function \l/ that leads to metal- 
lic or semiconducting behaviour depending on whether 
the transverse component of wavevector q is allowed to 
be zero [16, 17]. 



BOUNDARY CONDITIONS FOR THE 
EFFECTIVE MASS MODEL 

In order to obtain hard wall boundary conditions for 
the Dirac equation, we place an additional confinement 
potential at the boundary r = r^, 



—iva.V + cvMS (r — re) 



(10) 



of the vectors ni, 112, v\ and v-2,. In addition to p z we 
consider parity P = U x <g> I a , corresponding to a rota- 
tion by 7r about the z axis (a; — ► —x and y — > —y), and 
charge conjugation (C) and time reversal symmetry (T) 
that involve the complex conjugation operator combined 
with C = —iUy <8> <7 y and T = In ® cr y , respectively. 
The angles £ and £ mix terms with different symmetry 
with respect to p 2 : ( — and £ = correspond to even- 
ness with respect to p z whereas ( = tt/2 and £ = n/2 
correspond to oddncss. Since spin and/or helicity label 
different states at the same energy, values of £ and £ not 
equal to multiples of n/2 will lead to broken degeneracy. 
The angles A,T,0 and fl mix different symmetries with 
respect to combinations of P, C and p z . 



where c is a real constant and M is an arbitrary 4x4 
Hermitian, unitary matrix, M 2 = 1. The orientation of 
the boundary is defined by a unit vector normal to 
it, and we assume that the wavefunction is zero outside 
the confined region, but non-zero inside it. Then, wc 
integrate across an infinitesimal width of the boundary, 
giving 



(11) 



Substituting this equation back into itself, we find the 
requirements that c 2 M 2 = 1 (thus we set c = 1) and 

|nB.a, Af| = 0. The boundary condition can be writ- 
ten as * = where M = in B .aM, M 2 = 1 and 
{ne-Q, M} — 0, giving the result quoted in the introduc- 
tion Eq.(2). 

If, in the graphite coordinate system, we define 
the normal to the boundary as ne = (sin 77, cos 77, 0), 
where 77 is the chiral angle of the tube, then 
we may choose two mutually orthogonal 3D vectors 
as ni = (cos 77 sin £, — sin 77 sin £, cos £) and n 2 = 
(cos 77 cos £, — sin 77 cos £, — sin £) , and two additional or- 
thogonal 2D vectors as v\ — (cos£, sin£, 0) and z/2 = 
(— sin£, cos£, 0). This introduces two new mixing an- 
gles, £ and £: the arbitrary parameters contained within 
the boundary conditions describe the amount of mix- 
ing between different discrete symmetries. First, we 
note that the pseudo-spin of a 2D graphite sheet docs 
not transform in the same way as the spin of relativis- 
tic fermions, because certain transformations result in 
a swapping of the orientation of A and B atoms. This 
additional operation is described by a "pseudo-spin-flip" 
operator p z — W x <g> ia z that corresponds to a reflection 
in the (x, y) plane of relativistic fermions. For example, 
an active rotation of the 2D graphite sheet anticlockwise 
by 7r/3 about the perpendicular z axis, W (r') = (r), 
is described by Co = p z R(ir/3) = IT,, exp ((2wi/3)a z ) 
where R(9) = /nOexp ((i9/2)a z ) is a continuous rotation 
operator. 

Table 1 shows a summary of the discrete symmetries 
of the boundary conditions in terms of the orientation 



M 
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I n <8> ni.er 


n l = n (x,y) 


C — - 


-1 


+1 


+1 


-1 




ni = n z 


Cr. = 


+1 


+1 


+1 


-1 


U z ® n 2 .cr 


n 2 = n (x,y) 


(u =0 


+1 


-1 


-1 


-1 




n 2 = n z 


C — TL 
Stt 2 


-1 


-1 


-1 


-1 


u x .U®I a 


V\= I 


£« = 


+1 


+1 


+1 


+ 1 




v\= J 


t — TL 

Sit — 2 


-1 


-1 


+1 


-1 


V-2.YI (g> riB-Cr 


Vl= I 


£ — _2E 

Cm — 2 


-1 


+1 


-1 


-1 




vi= 3 


£« = 


+1 


-1 


-1 


+ 1 



Discrete symmetries of the boundary conditions 



SINGLE PARTICLE ENERGY SPECTRUM 

In order to illustrate the meaning of the general bound- 
ary conditions, we calculate the form of non-interacting 
single particle standing waves created by the boundary 
conditions and the corresponding energy spectrum. For 
simplicity, we will consider only metallic nanotubes with 
arbitrary chiral angle 77. We suppose that the x axis is 
perpendicular to the tube axis and we consider only the 
zero momentum transverse mode so that \E\ < 2-Kv/\C\ i \ 
where |Ch| is the circumference. The Dirac equation 
is diagonal in K-point space, so that, in the absence of 
boundary conditions, there are two right moving {^jP 



and ') and two left moving K and 
wave solutions: 
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where A, B, C and D are arbitrary constants, q is 
the wavevector along the tube and we consider q > 



(R) 

K 



and 



and E = svq, s = ±1. The solutions \I> 

are eigenvectors of pseudo-spin component S.J = 
in <S> e l1] ° z l' 1 <J v e~ v ' ]CTz l' 1 with eigenvalue +s, whereas the 
solutions and "f^' have eigenvalue — s. Also, 

the solutions and \I>^ are eigenvectors of pseudo- 
helicity -iS.V/|q| = H^In® e ir i a *l 2 {-ia y dy)e~ ir i a */ 2 
with eigenvalue +s, whereas the solutions and 



have eigenvalue — s. In the following we consider 
each of the four linear combinations Mi to M4 sepa- 
rately, and we consider a system with the same type 
of boundary condition on the right (at y = +L/2) 
and on the left (at y = —L/2). Wc introduce an 
index u = {R,L} = ±1 to label the right and left 
hand side so that the normal to the boundary, de- 
fined with respect to the graphite coordinate system, 
is iib = w(sin77,cos?7, 0), and we take into account the 
possibility of different mixing angles, A U ,T U ,Q U and 
and vectors ni = (ucos77sin£„, — usin^sin^cosCu), 
n 2 = (u cos 77 cos £„ , — u sin 77 cos , — sin£ u ), v\ — 
(cos£ u ,sin£ M ,0) and i/ 2 = (- sin£ u , cos£ u , 0). 



Ml: diagonal boundary conditions 

With the above definitions of the mixing angles, the 
boundary condition ^ = M\fy produces the following 
relations between the components of the wavefunction at 
the interface: 

u sin (( u + A u ) e~' in ip AK - [1 + cos (( u + A„)] ip BK = 0, 
u sin (C„ - A u ) e +tn ip AR - [1 - cos (C„ - A u )] %l> Bk 



0. 



The equations are diagonal in K-point space so do not de- 
scribe intervalley scattering. With these boundary con- 
ditions on the right (at y = +L/2) and on the left (at 
y = —L/2), standing waves at K are created from com- 
bining and and are labelled by helicity A = +s, 
and those at K are created from vj/^' and vfr^ and have 



label A 



-s. We find that 



K 



' K 



B = {-l) Pl Aexp [is {( R -( L )/2 + is (A R -A L )/2], 
D = {-If 2 C exp[-is (C h -Cl)/2 + is [A R -A L )/2], 

and the corresponding wavevectors are 

(A=+ S ) s(( R + ( L ) _ s(A R + A L ) Trpi 

q ' 2L 2L L ' 

_(A=-a) 1 s (Cr + Cl) s(A R + A L ) W P2 

q ' + 2L 2L L ' 

where {pi,P2} are integers such that q > 0. Using E = 
svq shows that the mixing angles £r and (l break K- 
point degeneracy whereas A R and A^ break electron-hole 
symmetry. 



M2: off-diagonal boundary conditions 

The boundary condition = M 2 V I' is equivalent to 
the following relations between the components of the 
enveloped wavefunction at the interface: 



ipAK + u sin T„e+ M '- J «" V„ k - cos T u e- l «" 



0. 



tpBK - MsinT u e %ri * €u V>b.k ~ cos T « e ^^ak = °- 

The equations are off-diagonal in K space so describe 
intervalley scattering. Standing waves are created from 
combining and with spin eigenvalue £ = +s, 

and and with spin eigenvalue £ — — s. We 

find that 

D = (-If 1 A C xp[ iS (T R - T L )/2 + is (& + &)/2], 
B = (-l) P2 Cexp[is(r R -r L )/2-is(Z R + Z L )/2]. 

and the corresponding wavevectors are 

„(£=+ s ) = s(T R + T L ) (Zr-Zl) | *Pi 



2L 



2L 



_(E=-«) s(T R + T L ) t , irp 2 

q ' 2L 2L L ' 

where {pi,P2} are integers such that q > 0. The angles 
£r and £l break degeneracy whereas Tr and Tl break 
electron-hole symmetry. 



M3: mixed boundary conditions (i) 

The boundary condition ^ = produces the fol- 

lowing relations between the components of the wave- 
function at the interface: 

^ak (1 - sin n u cos ( u ) + u cos O u e +l?7 ~ ,€ " tp A ^ 

-usinf2„ sin (; u e +%r, ipBK = 0, 
ipBK (1 + sinQ u cosCu) - ucosl]„e^' rl ^^ 

— usmfl u sm( u e~ lri ipAK = 0. 

The matrix M 3 has both diagonal and off-diagonal in K- 
point space parts. Standing waves are created from linear 
combinations of all *^ } , ^ { ? ] and We find 

that 

B = sin fi M e l;s " c " Ae iuqL + isucosflue'^ Ce mqL , 
D = sintt u e- lsu ^Ce luqL + isu cos fl u e +t ^ Ae tuqL . 

and the corresponding wavevectors are given by 

cos (2qL) = cos/3; 
cos/3 = sinfi^sinfii cos (( R + (l) 
- COS tl R COS cos - ^l) , 
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where q > and < (3 < tt. The energy levels are 

where n = {1,2,3,...} are integers. The spectrum al- 
ways has positive-negative energy symmetry, but broken 
degeneracy for j3 ^ {0, n}. 

M4: mixed boundary conditions (ii) 

The boundary condition ^ = M 4 \I> produces the fol- 
lowing relations between the components of the wave- 
function at the interface: 

ipAK (1 + sin 6„ sin ( u ) - u sin 9„ cos C u e +%r, ipBK 

— cosQ u e u 4> B K = 0, 
ipBK (1 - sin Q u sin ( u ) - u sin 9„ cos ( u e~ lv ipAK 

-<x>sQ u e~ % ^ip A ^ = 0. 

The matrix M4 has both diagonal and off-diagonal in 
K space parts. Standing waves are created from linear 
combinations of all * ( *\ * ( ~ L) , ¥~ R) and tf^ 5 . Wc find 

K ' K ' K K 

that 

B = isusixi<d u e t8U ^Ae luqL +cos<d u e- l ^Ce mqL , 
D = isusme u e- isu ^Ce iuqL + cosQ u e +l ^ Ae luqL . 

and the corresponding wavevectors are given by 

cos (2qL) = cosk; 
COSK = COS 9#COS Ql cos (£r — 
- sin Si? sin 6l cos (Cr + (l) , 

where q > and < k < ir. The energy levels are 

s ^^{ 2 T'- ± 2t}' (13) 

where n = {1,2,3,...} are integers. The spectrum al- 
ways has positive-negative energy symmetry, but broken 
degeneracy for n ^ {0, 7r}. 

DISCUSSION 

In this paper, we considered the effective mass model 
of spinless electrons in single wall carbon nanotubes that 
describes slowly varying spatial envelope wavefunctions 
$ (q, r) with small wavevectors q in the region of lin- 
ear dispersion E = ±w|q| near the K-points. Taking 
into account the two inequivalcnt K-points, the envelope 
wavefunctions "3/ obey the Dirac equation for masslcss 
fermions, written in terms of four component spinors, 
with the role of spin assumed by the relative ampli- 
tude of the wave function on the sublattice atoms ( "A" 



and "B"). Wc found that energy independent hard wall 
boundary conditions for the Dirac equation may be writ- 
ten as ^ = where M is an Hcrmitian, unitary 4x4 
matrix M 2 = 1 with the additional constraint that it an- 
ticommutes with the component of the current operator 
normal to the boundary. All possible linear combinations 
of matrices M obeying these constraints were expressed 
in terms of a small number of arbitrary parameters, mix- 
ing angles, that describe mixing between boundary con- 
ditions with different discrete symmetries. Then, in order 
to illustrate how the presence of non-zero mixing angles 
breaks K-point degeneracy and electron-hole symmetry, 
we evaluated the resulting energy level spectra for non- 
interacting electrons in finite-length metallic nanotubes 
and the corresponding standing wave envelope functions. 

The intention of this paper was to classify all possible 
boundary conditions for the spatially long-range envelope 
functions of a closed nanotube with length much greater 
than its circumference L 3> L c . The analysis was re- 
stricted to energy independent boundary conditions, al- 
though in principle they could be generalised by perform- 
ing a gradient expansion. We focused on armchair tubes 
and did not consider the possibility of edge states that 
may exist in zigzag graphite edges [21, 22, 23]. Rather 
than modelling the microscopic details of a boundary, 
such as shape or roughness, the nature of the bound- 
ary is characterised by mixing angles that describe the 
degree of symmetry breaking. In practice, the correct 
choice of a particular set of boundary conditions and val- 
ues of symmetry mixing angles needed to describe a given 
nanotube will depend on experimental details and may 
not be known beforehand. To illustrate this, wc com- 
pare our results to the idealised microscopic models of a 
a capped armchair nanotube and a model of a bound- 
ary obtained by setting the wavefunction to zero along 
a straight line of atoms. We find that the description of 
a capped nanotube considered in Ref. [13] corresponds 
to our off-diagonal boundary conditions M2, 'f = M 2 ^, 
with T u = or ir so that the component ^ak of the wave- 
function is related to iPbk a ^ ^ ne boundary. A model of 
a boundary obtained by setting the wavefunction to zero, 
which is equivalent to a particle- in-a-box model [14], also 
corresponds to our off-diagonal boundary conditions M2, 
1 3> = M 2 ^, but with T u = ±tt/2 so that the compo- 
nent ipAK of the wavefunction is related to ip A j^ at the 
boundary. Although the two models are described by a 
different mixing angle T, they both correspond to the 
off-diagonal boundary conditions, introduce inter-valley 
scattering at the boundary and, in general, they break 
K-point degeneracy with the mixing angle £ u dependent 
on the length of the nanotube. 

The authors thank J T Chalker and C J Lambert for 
discussions, and EPSRC for financial support. 
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APPENDIX 

In this appendix, we briefly describe the method of 
finding linear combinations of matrices that satisfy the 
constraints on M described by the boundary conditions, 
Eq. (2). Any 4x4 matrix may be written as a linear 
combination of the matrices In ® I„, In ® (n.cr), (V.II) ® 
I ai and (y.n)<8> (n.cr), where n and v are arbitrary three- 
dimensional space-like vectors. The first step is to find 
all the linear combinations that produce the unit matrix 
In <8> 1(7 when squared: 

M a = In® h, 

Mb = cos 9 (In ® ni .a) + sin 6 (i>\ .11 ® n 2 .a) , 
M c = cos (^i.II ® I CT ) + sine/) (^2-11 (g) n 2 .cr) , 

where the vectors are unit vectors with additional con- 
straints ni.n 2 = 0, ^i.^2 = that ensure no cross-terms 
survive. The next step is to find the conditions under 
which the matrices M a , Mb, and M c anticommute with 
the operator riB.Q that is proportional to the component 
of the current operator normal to the interface. Clearly, 
M a does not anticommute, so it is discarded. The matrix 
Mb anticommutes if ni.ne = 0, n 2 .riB = 0, and v\ = k 
(Ml), or if ni.ne = 0, n 2 = ne, and v\ is confined to 
the (x,y) plane (M3). The matrix M c anticommutes if 
v\ and V2 are confined to the (a;, y) plane and n 2 = ne 
(M2), or if v\ is confined to the (x, y) plane, V2 = k, and 
n 2 .n B = (M4). 
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